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Dynamical fat link fermions 
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The use of APE smearing or other blocking techniques in fermion actions can provide many advantages. There 
are many variants of these fat link actions in lattice QCD currently, such as FLIC fermions. Frequently, fat link 
actions make use of the APE blocking technique in combination with a projection of the blocked links back into 
the special unitary group. This reunitarisation is often performed using an iterative maximisation of a gauge 
invariant measure. This technique is not differentiable with respect to the gauge field and thus prevents the use of 
standard Hybrid Monte Carlo simulation algorithms. The use of an alternative projection technique circumvents 
this difficulty and allows the simulation of dynamical fat link fermions with standard HMC and its variants. 



1. INTRODUCTION 

The use of smeared or "fat" links in lattice 
fermion actions has been of interest for some 
time now£Q. Fat Link Irrelevant Clover (FLIC) 
fermions have shown a number of promising ad- 
vantages over standard actions, including im- 
proved convergence properties [2] and O(a) im- 
proved scaling without the need for nonpertur- 
bative tuning 3 . Furthermore, a reduced ex- 
ceptional configuration problem has allowed ef- 
ficient access to the light quark mass regime in 
the quenched approximation [3]. However, recent 
progress in the field has shown that the behaviour 
of quenched QCD can differ from the true theory 
both qualitatively and quantitatively in the chi- 
ral regime |5I6| . Thus, interest is now focusing 
on dynamical QCD, be it (truly) unquenched, or 
partially quenched. 

In particular, as it is in the chiral regime where 
the difference from the valence approximation will 
be highlighted, we would like to go to light quark 
masses in dynamical QCD. This is an extremely 
computationally expensive endeavour. One might 
hope that the excellent behaviour at light quark 
mass displayed by FLIC fermions will carry over 
from the quenched theory to the unquenched one. 
This brings us to the issue of generating dynam- 
ical gauge field configurations with the fermionic 
determinant being that of the FLIC action. 



2. MaxReTr PROJECTION 

The standard technique for simulating dynam- 
ical fermions has for some time now been Hy- 
brid Monte Carlo (HMC) 7 . The HMC algo- 
rithm alternates between a global Metropolis ac- 
cept/reject step and a gauge field update through 
the use of a hybrid Molecular Dynamics integra- 
tion, 

U^x,t+At) = U f ,(x,T)exp(iATP ll (x,T)), (1) 



Pp(x,t+At) =P ix {x,t)-U, 1 {x,t) 



SS 



SUfj,(x,' 



Here, the Pfi(x) are the momenta canonically con- 
jugate to the gauge field. The update of the mo- 
menta involves the derivative of the action with 
respect to the gauge field, a quantity which is also 
made use of in efficient local update algorithms. 

When attempting to use HMC with FLIC 
fermions, or some other fat link action for that 
matter, it is the calculation of which has 
proven problematic. FLIC fermions are clover- 
improved fermions where the irrelevant operators 
are constructed using APE smeared links |8I9| . 
These fat links are constucted by performing sev- 
eral sweeps of APE smearing across the lattice. 
Each sweep consists of an APE blocking step, 



vW(x)[U^] = (l-a)^+^ri- 



U 



1 



2 



Sweep 



Unit Circle 



MaxReTr 



0.866138301214314 

1 0.960313394813806 

2 0.980735000838119 

3 0.988384926461589 

4 0.992103013943516 

5 0.994182852413813 

6 0.995457365275018 

7 0.996293668622924 

8 0.996878305318083 



0.866138301214314 

0.960348747275940 
0.980751346847750 
0.988393707639555 
0.992107844842705 
0.994185532052157 
0.995458835653863 
0.996294454083006 
0.996878710433084 



Table 1 

The mean link uq for a single configuration as a function of number of APE smearing sweeps, for the two 
different projection methods. The boldface indicates significant digits which match. The configuration is 
a dynamical gauge field with DBW2 glue and FLIC sea fermions, at f3 = 8.0, k = 0.128. 



(3) 



followed by a projection back into SU(3) of the 
links, uj?\x) = V(V^ n \x)). Typically, this pro- 
jection is performed by iteratively maximising the 
following gauge invariant measure, 

Ui n) (x)[U {n - 1) }= max ReTr(f/'(x)U (n)t (.T)). 



(4) 



which is special unitary 1 . This technique pro- 
vides an alternative method of SU(3) projection 
(breaking the Z3 ambiguity by choosing the prin- 
cipal value of the cube root). Using the mean link 
as a measure of the smoothness of the smeared 
gauge field, Table ^ indicates that the two meth- 
ods presented here are nearly equivalent. 

However, the matrix inverse square root func- 
tion can be approximated by a rational polyno- 
mial (whose poles lie on the imaginary axis), 



Any fermion action which makes use of this re- 
unitarisation method cannot be simulated using 
HMC, as the use of the iterative algorithm pre- 
vents one from being able to construct 

3. UNIT CIRCLE PROJECTION 

Given any matrix A, then X'X is hermitian 
and may be diagonalised. Then it is possible (for 
det X 7^ 0) to define a matrix 



W 



X 



vxxx 



(5) 



whose spectrum lies on the complex unit circle 
and is hence unitary. Furthermore, W possesses 
the same gauge transformation properties as X. 
We can then constuct another matrix, 



w = 



1 



VdetW 



-.W 



(G) 



W[X] a W k [X] = d X(X^X+c )J2 



1=1 



At A - 



- Q 
(7) 



This approximation is differentiablc in a matrix 
sense for all A for which the inverse square root 
can be defined. This means that we can construct 
J|j for fermion actions which involve unit circle 
projection. 

4. RESULTS 

Having now written down the APE smearing 
prescription (with projection) in a differentiable 
closed form, the equations of motion necessary 
for the use of the HMC algorithm can be de- 
rived. Even in the simple case of APE smearing 
these equations are quite complex, but we will say 
something here about the mechanism for deriving 

lr The necessary cube root is absent in Ref. 1101 
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them. Given the FLIC fermion action <Sflic(^0j 
which has explicit dependence upon the thin links 
U and the fat links U', it is straightforward to 
write down the partial matrix derivatives and 
■Sjjr- However, to write down in one step the to- 
tal matrix derivative ^ would be difficult, and 
inefficent numerically. Instead, one can use some 
calculus to both simplify the process and make 
it numerically efficient. This is done through the 
use of the matrix chain rule, 

dS_ _ dS_ dS_ dlT_ 

dU ~dU + dU>~ k ~dU' ( ' 

The full details of the derivation are presented 
elsewhere |11|. We find that the numerical over- 
head for the equations of motion are small, and 
the generation of dynamical gauge fields with 
FLIC fcrmions is still dominated by the neces- 
sary conjugate gradient inversion required as part 
of the pseudo-fermion formulation. Results are 
given in Table [3 

5. CONCLUSION 

Unit circle projection can be written in closed 
form and represented by a rational polynomial 
approximation. The matrix rational polynomial 
is differentiable with respect to the gauge field 
and thus enables one to derive the equations of 
motions necessary for the use of HMC with FLIC 
fermions 2 . This provides a significant advantage 
over the MaxReTr method, where the absence of 
a differentiable form precludes the use of HMC 
and hence the availability of a simulation algo- 
rithm that scales linearly with the lattice volume 
V, although there are 0(V 2 ) alternatives |13| . 
This means that we now have an 0(V) algorithm 
available for the dynamical simulation of FLIC 
fermions. Furthermore, the method is general 
and can be applied to any fermion action with 
reuniterisation, including overlap fermions with 
a fat link kernel |14I15I16I17| . or other types of 
fat link actions ^2] that may involve alternative 
smearing techniques |T5] . 
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Table 2 

The step sizes, acceptance rate, mean link, lattice spacing and pion mass (in GeV) for different dynamical 
simulations, using (two flavours of) FLIC sea quarks and both Lusher- Wiesz (IMP) glue and DBW2 glue. 
These results are obtained from 20 12 3 x 24 configurations. Simulations are done using multiple time 
step HMC with trajectories of unit length. The lattice spacing is set by ro via the static quark potential. 



